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FOREWORD 

This is the final report for Phase II for Contract NASI- 18465 
Processed Synthetic Aperture Radar (SAR) Data, sponsored by NASA Langley 
Research Center (LaRC). This report is second in a series of 
statistical reports. The overall effort has as its thrust the 
description of the ground clutter at an airport and in the surrounding 
area. In Phase I of this activity, SAR data of airports which existed 
in the ERIM SAR data archive were examined for utility to this program. 
Eight digital scattering coefficient images at high resolution and 
coarse resolution were created. The coarse resolution images were 
provided to NASA LaRC for use in their Microburst/Clutter/Radar 
simulation program and the high resolution images underwent a 
statistical clutter analysis at ERIM. In this phase of the program, SAR 
data were collected on an opportunity basis at the Philadelphia Airport 
using a set of radar parameters which more closely matched those which 
are anticipated to be encountered by an aircraft on its approach to an 
airport. This report describes these data and the results of the 
clutter study. During Phase III, a dedicated SAR mission was flown of 
the Denver Stapleton Airport and surrounding area. A wide variety of 
geometries and scene contents were acquired. These data and study 
results are presented in the third report. 

The work reported here was performed by members of the Radar Science 
Laboratory, Advanced Concepts Division, Environmental Research 
Institute, under the direction of Dr. S.R. Robinson. The principal 
investigator for this project was Dr. R.G. Onstott. The contract was 
monitored by E.M. Bracalente, NASA Langley Research Center, Hampton, 
Virginia. 

i i i 


preceding page blank not filmed 


Wis e /( I NTEfmQHAhtl BLANK 


The authors gratefully acknowledge the support of the Radar Science 
Laboratory Staff during the project, in particular, Mr. Christopher 
Wackerman, Mr. Steven Lerman, Ms. Rita Woods, Ms. Janice Anquetil, Ms. 
Marsha Allen, and Mr. James Schipp. 



TABLE OF CONTENTS 


FOREWORD iii 

LIST OF TABLES vi 

LIST OF FIGURES vi i 

I. INTRODUCTION 1 

II. DATA COLLECTION 5 

III. IMAGE DESCRIPTION 9 

IV. PROCESSING AND CALIBRATION 11 

V. CLUTTER ANALYSIS PROCEDURES 13 

VI. CLUTTER ANALYSIS RESULTS 15 

VII. CONCLUSION 23 

References 81 

Appendix A A-l 

Appendix B B-l 

Appendix C C-l 

Appendix D D-l 


v 


LIST OF TABLES 


1. Comparison of Preferred and Actual Radar and Geometry 

Parameters 24 

2. Image Composition Areal Analysis 25 

3. Statistical Summary of Backscatter Cross-Sections 26 

4. Hard Targets Represented as Sigma Values 27 


LIST OF FIGURES 

1. SAR Imaging Geometry 28 

2. SAR Imaging Geometry 29 

3. P-3 SAR Swath Modes 30 

4. SAR Data Collected of Philadelphia Airport and Surrounding 

Area on 7 October 1988 3 ^ 

5. Bar Chart Presentation of Scattering Coefficients 

Values in the Philadelphia Airport Image 32 

6. Histogram of Radar Scattering Coefficients for the 

Philadelphia Image 33 

7. -30 dB Threshold 3 ^ 

8. -20 dB Threshold 35 

9. -10 dB Threshold 36 

10. 0 dB Threshold 37 

11. 10 dB Threshold 38 

12. Distribution of Radar Scattering Cross Section for the Whole 

Airport Image 3 ^ 

13. Residential Clutter Distribution (40° to 49° incidence angle)... 40 

14. Residential Clutter Distribution (50° to 59° incidence angle)... 41 

15. Residential Clutter Distribution (60° to 64° incidence angle)... 42 

16. Residential Clutter Distribution (65° to 69° incidence angle)... 43 

17. Residential Clutter Distribution (70° to 74° incidence angle)... 44 

18. Forest Clutter Distribution (40° to 49° incidence angle) 45 

19. Forest Clutter Distribution (50° to 59° incidence angle) 46 

20. Forest Clutter Distribution (60° to 64° incidence angle) 47 

21. Forest Clutter Distribution (65° to 69° incidence angle) 48 

22. Forest Clutter Distribution (70° to 74° incidence angle) 49 


LIST OF FIGURES (cont.) 

23. Forest Clutter Distribution (75° to 79° incidence angle) 50 

24. Grass Clutter distribution (40° to 49° incidence angle) 51 

25. Grass Clutter Distribution (50° to 59° incidence angle) 52 

26. Grass Clutter Distribution (60° to 64° incidence angle) 53 

27. Grass Clutter Distribution (65° to 69° incidence angle) 54 

28. Grass Clutter Distribution (70° to 74 ° incidence angle) 55 

29. Grass Clutter Distribution (75° to 79° incidence angle) 56 

30. Grass Clutter Distribution (80° to 84° incidence angle) 57 

31. Grass Clutter Distribution Including Values at Image 

Minimum (75° to 79° incidence angle) 58 

32. Grass Clutter Distribution Including Values at Image 

Minimum (80° to 84° incidence angle) 59 

33. Water Clutter Distribution (40° to 49° incidence angle) 60 

34. Water Clutter Distribution (50° to 59° incidence angle) 61 

35. Water Clutter Distribution (65° to 69° incidence angle) 62 

36. Water Clutter Distribution (70° to 74° incidence angle) 63 

37. Water Clutter Distribution (75° to 79° incidence angle) 64 

38. Water Clutter Distribution (80° to 84° incidence angle) 65 

39. Incidence Angle vs. Scattering Coefficient Plot for 

Forest Clutter 

40. Incidence Angle vs. Scattering Coefficient Plot for 

Residential Clutter 67 

41. Incidence Angle vs. Scattering Coefficient Plot for 

Grass Clutter gg 

42. Incidence Angle vs. Scattering Coefficient Plot for 

Water Clutter 69 

43. Polynomial Fit to Forest Data 

vi i i 



LIST OF FIGURES (cont.) 


44. Polynomial Fit to Residential Data ,L 

45. Polynomial Fit to Grass Data 72 

46. Distribution Plot for the Airport Terminals 73 

47. Distribution Plot for Building Clutter 74 

48. Distribution Plot for Refinery Tank Clutter 75 

49. Distribution Plot for Dock 76 

50. Distribution Plot for Bridge Clutter 77 

51. Distribution Plot for Plane Clutter 78 

52. Distribution Plot for Barge Clutter 79 

53. Distribution Plot for Parking Lot Clutter 80 




I. INTRODUCTION 


Low altitude microburst windshear represents a significant hazard to 
aircraft, particularly during take-off and landing; the intense down 
drafts and the resultant divergent outflow, can have a significant 
effect on the characteristics of the endangered aircraft. When 
encountered at low altitude, the pilot has little time to react 
correctly to maintain safe flight. The Federal Aviation Administration 
(FAA) , jointly with the National Aeronautics and Space Administration 
(NASA), has sponsored an investigation into the development of airborne 
sensors to detect microburst windshear. One sensor of interest is the 
microwave Doppler radar operating at X-band or higher frequencies. 
Critical to the analysis of the capability of such a sensor to perform 
this detection is the microwave backscatter description of both the 
microburst event and the clutter background, especially during the 
approach and departure from an airport. The Environmental Research 
Institute of Michigan (ERIM) has provided NASA Langley Research Center 
(LaRC) with calibrated high-resolution synthetic aperture radar (SAR) 
images of selected airport scenes from the ERIM SAR data archive for use 
in their Microburst/Clutter/Radar simulation program and to 
statistically characterize the ground clutter surrounding these 
airports. In addition, statistical analyses of these airport 
environments have been performed by ERIM to describe the range of 
scattering conditions encountered. 

NASA LaRC has developed this Microburst/Clutter/Radar simulation 
program to assess the performance of Doppler radars for this application 
and to test potential signal processing techniques. The simulation 
program models the output of an airborne Doppler radar as it views a 
low-level microburst along the approach to an airport. Inputs to this 
simulation include airport ground clutter data base, a simulated 
microburst data base, the operating parameters of the proposed weather 
radar, and candidate signal processing software for use in detection. 

In the operation of the simulator the received signal amplitude level 
for each range bin is calculated. Each range bin include contributions 



from both the microburst and the ground clutter. 

Many questions arose during this study. What is the general 
description of the airport clutter environment? How does this 
environment change from airport to airport? How complex does an airport 
scene have to become before it degrades detection? In general, the 
limits of detectability of microburst events must be established for an 
airborne Doppler radar to be an effective tool in the prevention of 
aircraft windshear catastrophes. ERIM has performed a clutter analysis 
of the Philadelphia Airport to extend the work done in Phase I to 
provide significant information for establishing limits of 
detectability. Clutter types, mean backscatter intensities, probability 
distributions, and areal extent of the clutter types in the image have 
been determined. With the analysis of the statistical characteristics 
of clutter in the scene, the effects of scene composition were studied. 
These results may then be extrapolated to describe the clutter at 
candidate airports. 

The analysis of the eight images from ERIM's archives in Phase I of 
this three phase program has provided important information concerning 
clutter returns from airport scenes. Clutter backscatter responses for 
the same clutter type were found very similar. Differences in their 
means were attributable to differences in incidence angle. Probability 
density functions which describe the scattering of a particular clutter 
type were nearly identical from image to image. This is an important 
result because it indicates a high degree of stability in the returns 
which would be expected from clutter around an airport. A Doppler radar 
designed to detect Microburst windshear at one particular airport may 
not have to compensate radically for different clutter scenes as the 
plane travels from airport to airport. Through the analysis of these 
eight images we gained information about the types and characteristics 
of radar clutter surrounding airports. However, further analysis of 
strongly backscattering clutter regions was needed. The data analyzed 
in Phase I indicated that although only a fraction of the clutter 
surrounding the airport is of an intensity to be of concern, almost all 
of it is located on the airport grounds. It was concluded that detailed 
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analysis of the clutter from specific types of buildings, planes, and 
other airport vehicles should be performed, and that the motion of 
vehicles and planes needs to be examined. 

As anticipated, the data in the SAR archive, by itself, did not have 
the breadth to fully describe the range of airport clutter scenes 
possible. This, of course, was the impetus behind Phases II and III, 
which was to collect additional clutter data on a not-to-interfere basis 
during the calibration flights of a new SAR sensor, and to conduct a 
dedicated SAR data collection. The SAR archive allowed us to begin a 
statistical survey of the clutter environment at a few selected 
airports. This information was enriched in Phase II by gathering data 
at smaller depression angles. It was also determined by NASA LaRC 
personnel that an airport scene of 13 km x 13 km was optimum. In Phase 
I images were limited to sizes of about 6 km x 13 km. These scenes were 
then filled in by a variety of techniques employed by NASA LaRC 
personnel. During the Philadelphia collection usable data over a 13 km 
x 13 km region were acquired. 

During this phase the goal was to make use of the opportunity to 
acquire airport clutter during the test and evaluation flights of a new 
radar, commonly referred as the NADC/ERIM P-3 SAR, designed for remote 
sensing applications. The goal was to collect data on a not-to- 
interfere basis using the geometry and radar parameters specific to the 
Wind Shear Detection clutter problem. A number of opportunities 
presented themselves, but did not result in the desired collection of 
data (i.e., because of weather, system problems, and tight schedules). 
Data were collected of the Key West Naval Air Station (15 January 1988), 
Willow Run Airport (17 May 1988), and Whidbey Island Naval Air Station 
(16 July 1988). The maximum incidence angle for these images was 74° 
and the clutter was not extensively urban so these data were not 
exploited. On 7 October 1988, SAR data of the Philadelphia Airport 
which is located in a densely populated area were collected using the 
geometry and primary radar parameters that were judged the most optimum. 
The analysis of these data is the object of the discussion of this 
report. 
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II. DATA COLLECTION 


The radar used during this data collection operates at frequencies 
of 1.25, 5.26, and 9.375 GHz. It is fully polarimetric, able to record 
4 complex images, operates in the strip map mode, and is capable of 
producing 1.6 meter range by 2.2 meter azimuth resolution images. Four 
high density recording channels are available, each capable of recording 
4096 range elements. The operator is able to select wavelength, 
transmit polarization, receive polarization, fine or coarse (1.6 to 3.2 
m) , range resolution and slant range for each channel. It is possible 
to operate in a wide swath mode so that a 32 km strip map is made at one 
wavelength and one polarization or the four channels can be used to 
record at one frequency, all four linear polarizations, and fine 
resolution over a 6 km range interval. 

A more detailed description of the SAR is provided in Appendix A. A 
data collection plan was developed to provide SAR data which closely 
matches the geometry of an aircraft on its landing approach. The 
various radar parameters and geometries are addressed individually 
below. The preferred and actual SAR configurations are detailed in Table 
1. The planned SAR imaging geometry is provided in Figures 1 and 2. 

Frequency 

The P-3 SAR operates at 1.25, 5.25, and 9.375 GHz. The frequency of 
9.375 GHz (X-band) was selected because it corresponds with existing X- 
band weather radars which have proven useful in detecting weather 
related features. 

Polarization 


The P-3 SAR Facility is fully polarimetric. This means that up to 
four complex images at VV, VH, HV, and HH polarizations may be recorded 
at the selected frequency of operation. The first letter indicates the 
transmit polarization and the second the received polarization. The 
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letter V is for vertical and H is for horizontal. With the recording of 
these four complex images additional images at arbitrary transmit and 
receive polarizations can be synthesized (i.e., circular polarization). 
Recording the four complex images also allows for the complete 
description of the scattering properties (i.e., by providing the cross- 
correlation information between the scattering coefficients) of the 
various elements within the clutter scene. Such information is critical 
to the understanding of the scattering mechanisms at play, and 
potentially to the determination as to whether troublesome scattering 
scenes may be suppressed. During this data collection, priority was 
given to recording four complex, 12 km x 12 km images (VV, VH, HV, and 
HH). These were produced in one pass using coarse resolution. 

Depression Angle 

An aircraft during its final approach is typically flying on a 3 
degree glide slope. To simulate this glide slope during the SAR data 
collection, it was considered optimum if images were centered about a 3 
degree depression angle. Depression angle is defined as the angle off 
the horizon and is equal to 90°-incidence angle. 

Resolution 


For any swath modes recorded, the P-3 SAR can operate either in high 
resolution, with an azimuth resolution of 2.8 m, and a range resolution 
of 1.6 m, or in low resolution, with an azimuth resolution of 2.8 m and 
a range resolution of 3.0 m. Using low resolution, the range coverage 
of each channel is 9830 m, and using high resolution, the range coverage 
of each channel is 4915 m. Since the goal here was to produce a 13 km x 
13 km image about the airport, coarse resolution was chosen so that 
four complex 12 km x 12 km images (VV, VH, HV, and HH) may be produced 
in one pass. Operation in this mode was selected since the final images 
would be digitally processed to a resolution of 20 m. 
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Image Size 


An image about 13 km x 13 km with the airport touchdown point near 
the center of the image was considered optimum by NASA LaRC personnel. 
The system may be operated in three different swath modes. The 
configuration of these swath modes is illustrated in Figure 3. In 
single-swath mode (Figure 3a), the radar collects the same range 
coverage in each channel using four different frequency and polarization 
combinations. In double-swath mode (Figure 3b), twice the range 
coverage of single-swath mode is recorded but only two frequency and 
polarization combinations can be used. In quad-swath mode (Figure 3c), 
four times the range coverage of single swath mode is recorded but only 
one frequency and polarization combination can be used. In double- and 
quad-swath mode there is no range overlap between channels that are 
contiguous in range. Images of 12 km x 12 km may be formed in 1 pass 
using the double swath-mode and was selected. 

Airport/SAR Orientation 

To simulate the geometry of what a doppler radar will see on final 
approach, the SAR needed to be oriented such that the path of the 
aircraft was perpendicular to the main runways of the airport. If 
additional passes were possible then alignment with secondary runways 
was proposed to provide complementary views of the airport from which 
aspect angle dependence would be studied and the composition (i.e., 
grass, forest, urban, etc.) of the airport scene may be supplemented. 

Airport Selection 

The P-3 SAR flights made during this period were strongly driven by 
the desire to fully document sensor performance. It was anticipated 
that there would be numerous opportunities to gather these data on an 
opportunity basis. It was originally planned to establish calibration 
arrays at Willow Run Airport which were to serve as the primary test 
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site for the SAR calibration program. During the planning stage of this 
phase it was not known that opportunities at other airports would be 
possible. The need to examine a variety of airports with their unique 
clutter environments was recognized. 
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III. IMAGE DESCRIPTION 


The SAR data shown in Figure 4 were collected of the Philadelphia 
International Airport on 7 October 1988. Using a frequency of 9.35 GHz 
and a horizontal transmit and receive polarization. The radar was flown 
at an altitude of 5287 feet and grazing angles subtended by the image 
ranged from 7.8° to 52.2°. The radar was flown with a heading of 180 
degrees and a look direction of 270 degrees. North is at the left of 
the image. 

The Philadelphia International Airport is located near the center in 
the far range of the image. The main terminal and concourses can be 
seen in some detail, and some of the runways are visible. To the south 
of the airport is the Delaware River. Ships were observed in transit 
and at the far east (top) portion of the image a ship wake is observed. 
South of the river is rural southern New Jersey. Most of this area 
appears to be parks, open fields, or residential areas. On the far 
south eastern shore of the Delaware River is an oil refinery. One of 
the arrays of bright returns is the refinery whereas the other arises 
from a field of storage tanks. The strong returns which extend into the 
river near the refinery are the shipping terminals. 

Just north of the airport is Interstate 95. Following the 
interstate east from the airport is Fort Mifflin U.S. Army Base, which 
lies at the junction of the Delaware and Schuylkill Rivers. The very 
bright area and the apparently grassy area east of the Schuylkill River 
and south of Interstate 95 are both part of the Philadelphia Naval 
Shipyard. 

Directly north of the interstate and the shipyard are some 
interesting urban features. The golf courses and lakes of F. D. 
Roosevelt Park appear as darker areas just north of the bright part of 
the shipyard. To the east of the park is a sports complex which 
consists of Veteran's Stadium, J. F. Kennedy Stadium, and Spectrum 
Sports Arena. Parking areas surround these facilities. Lightpoles 
within these lots can be seen. To the north of the sports complex and 
extending to the north is downtown Philadelphia. 
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IV. PROCESSING AND CALIBRATION 


Processing and calibration of the image proceeded in much the same 
manner as the images processed from the SAR archive, the object of the 
work done during Phase I. Refer to Volume I of this series of reports 
(Onstott and Gineris, 1988) and to Appendix B for a description of this 
calibration procedure. The procedure deviated slightly due to recent 
modifications in the processor which enabled us to produce an image with 
the full ground coverage desired, instead of an image produced from the 
mosaicking of multiple images. The resultant phase history tape used 
for the Philadelphia image represents a 8192 element by 4096 element 
file. With a pixel spacing of 1.62 m in azimuth and 2.4 m in range, the 
Philadelphia image is produced from a reflectivity map of a 13,271 m by 
9830 m area. 

As with the archival images, the phase history of the Philadelphia 
image was focused in azimuth and range. This was achieved by convolving 
the data with a match filter of the transmitted radar chirp in azimuth 
and range. The data were processed to remove the effects of system 
noise and were then radiometrical ly corrected to compensate for the 
effects of range fall-off, the antenna gain pattern, and resolution cell 
power. The radiometrical ly corrected image was then converted to 
normalized radar cross-sections (NRCS) by normalizing the magnitude of 
the scattering cross-sections by the resolution area. The terms 
scattering coefficients and NRCS will be used interchangeably. 

The absolute calibration of the Philadelphia Airport image was 
performed based on data collected at the Patuxent Naval Air Test Center, 
Maryland, two weeks prior to the Philadelphia collection. For a radar 
operating in its linear region, a linear relationship will exist between 
the measured intensity of a point target in an image and the expected 
value of the backscattering cross-section of the target. The slope of 
the function is unity and the y-intercept of the function is a measure 
of the system gain function. The system gain function determined from 
the Patuxent data was then applied to the Philadelphia collection by 
accounting for short term variations in receiver gain and transmitted 
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V. CLUTTER ANALYSIS PROCEDURES 


In establishing the criteria for which microburst events can be 
detected, a careful statistical characterization of airport clutter is 
required. The questions that arose immediately were: what are the types 
of clutter which commonly occur at an airport, how does the clutter 
environment change with incidence angle and polarization, and what is 
the density and location of strong returns at an airport? The 
Philadelphia data were examined to address these questions. This 
analysis was performed on 4096 element x 4096 element, slant range a° 
image with the finest resolution possible to allow precise sampling of 
both distributed and point targets. This image has one independent 
sample per resolution cell. 

The first step in our analysis was to characterize the area 
surrounding the Philadelphia International Airport. This area is a 
mixture of industrial, suburban, and rural communities. Two analyses 
were performed to quantify the amount and characteristic clutter return 
of each type of community in the image. The first analysis performed 
was the differentiation of the image into its various constituent 
communities. The results of this areal analysis were tabulated and the 
percentage of each community within the image calculated. A second 
analysis was performed in which the image was divided into a° clutter 
groups at intervals of 5 dB. This enabled us to locate and quantify 
sources of similar backscatter intensity. 

Using standard statistical techniques, we performed an analysis of 
each major clutter type as well as the unique clutter areas. Regions of 
critical clutter types were located and the means, standard deviations, 
and coefficients of variation were calculated using the techniques 
described in Appendix C. Coefficients of variation were determined to 
examine the uniformity of these subregions. Most distributed targets 
have low coefficients of variation indicating the backscatter intensity 
does not vary widely from the mean value. Man-made targets typically 
have high coefficients of variation. The few pixels of very high 
returns are often times embedded in a background of low returns. The 
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coefficient of variation also provides textural information about a 
distributed targets. For example, a field of grass (a very uniform 
distributed target) is expected to have a lower coefficient of variation 
than a forest because forest returns may arise from dissimilar scatters 
such as trunks, crowns and the soil surface. During Phase I, skewness, 
kurtosis and probability density functions were determined indicating 
that most clutter types were well described using a gamma density 
function. These calculations were not further duplicated in Phase II. 

Regions of similar clutter properties were then employed in general 
clutter characterization and in the examination of their response with 
incidence angle. Areas of similar clutter type and incidence angle were 
merged. Five degree incidence angle bins were used in the range of 60 
degrees and beyond, since backscatter returns from distributed targets 
generally changes rapidly in this region. Ten degree incidence angle 
bins were used in the range from 30 to 60 degrees, where backscatter 
values are relatively constant. Histograms, means, standard deviations, 
and coefficients of variation were calculated. It should be noted that 
in order to compare the expected scattering cross-sections from point 
and man-made targets to normalized radar cross-sections the area extent 
of the target must be taken into account. 

Incidence angle effects in the data were examined by plotting the 
mean return of each clutter sub-region as a function of the mean 
incidence angle and are important in the analysis of both distributed 
and man-made targets. 
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VI. CLUTTER ANALYSIS RESULTS 


The image composition areal analysis of the Philadelphia 
International Airport and surrounding area is presented in Table 2. 

About 33% of the image is composed of water, i.e., the Delaware and 
Schuylkill Rivers. Another 18% of the image is comprised of fields, 
forests, and croplands. The majority of these areas lie to the south of 
the Delaware River in New Jersey, although some of the 'field' areas 
considered in the analysis are part of the park/army base/naval shipyard 
complex directly to the east of the airport. Industrial complexes (east 
of the airport) make up approximately 14% of the image content and the 
airport grounds itself occupy about 6%. About 16% of the far range 
image has scattering coefficients which are less than about -35 dB. The 
industrial areas to the east of the airport will be shown to represent 
the greatest clutter hazard, especially since it lies almost directly in 
the path of the major runways at the airport. 

Table 3 and Figure 5 present the mean, minimum value, maximum value 
and standard deviation for the scattering coefficients of each of the 
major clutter types in the SAR image. The location and statistics of 
individual clutter areas are given in Appendix D. From the bar chart 
provided in Figure 5 we see that, with the exception of the parking lot, 
every man-made hard target clutter group has a NRCS of -8 dB or greater. 
This correlates well to the results seen previously in the archival 
data. Distributed targets NRCS values are -20 dB and below. The 
parking lot filled with cars has a mean of -15 dB and does not fit into 
either of these two categories. This may be due to the composite nature 
of the parking lot. The coefficient of variation for the parking lot 
clutter area has a value of 15, one of the largest coefficients of 
variation in all of the data sets. A mixture of very large NRCS values 
, as would be expected from vehicles, and very low, as would be expected 
from a concrete or asphalt surface, would produce a mean NRCS value 
somewhere between the two extremes and a large coefficient of variation. 
For this image the mean NRCS was about -15 dB, with the largest NRCS 
measured of about 34 dB. The coefficient of variation has a value of 
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66.0, indicative of the very diverse nature of this image. 

Given that the scattering coefficients of the image tended to 
separate into either man-made or natural distributed-target the 
thresholding of the image allowed us to further determine the location 
and percentage of specific types of clutter. A histogram of radar 
scattering coefficients for the entire image and for 5 dB bins is 
provided in Figure 6. The images shown in Figures 7 through 11 were 
produced by setting individual pixels whose values are below the stated 
threshold to zero (black areas in the images) while pixel values equal 
to or greater than the threshold are set to 255 (white areas in the 
images). It should be noted that this operation has been performed on 
the compressed ( = 19 m resolution) version of the calibrated image. 
Because image compression involves averaging the distribution of 
coefficients, the resulting range of values is smaller, as expected, and 
the lowest and highest resulting coefficients will differ from those 
presented elsewhere. 

Based upon the bar chart shown in Figure 6, we see that fifty-four 
percent of the image has a value of -28 dB or less. Figure 7 reveals 
that these low return areas arise primarily due to the water and 
unclassified areas of the image. Low returns also exist in the 
residential and parking lot areas. Thirty-four percent of the image has 
NRCS values between -18 db and -28 dB. These areas include almost all 
of the distributed targets. Figure 8 shows all returns which are 
greater than -20 dB. The natural targets which produce NRCS values 
greater than -20 dB are shorelines, treelines, and certain tree returns. 
All other returns appear to be manmade in origin. This behavior was 
also suggested in the bar chart shown in Figure 5. The -15 dB to -20 dB 
NRCS value range appears to demarcate natural distributed returns from 
hard-target returns. At the -10 dB threshold, Figure 9, most of the 
returns which are visible come from hard-target sources. Most prominent 
are the returns from the naval shipyard, the airport terminal, the oil 
refineries, boats, and docks. A very small percentage of the returns 
with values above -10 dB come from the 'residential' areas in New Jersey 
south of the naval shipyard. These returns may actually be coming from 
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homes and other buildings dispersed throughout this area. Only 1.2% of 
the image has values of 0 dB or above and the majority of these returns 
are localized at the naval shipyard, the airport, and a few developed 
locations along the river and south of the airport as shown in Figure 
10. These clutter scenes produced the strongest returns in this entire 
image and are directly in the line of approach of the major airport 
runways. This is further illustrated in Figure 11 where the image 
thresholded at the +10 dB level. 

Histograms of radar scattering coefficients were calculated for the 
various clutter scenes as well as for the image as a whole and are 
presented in Figures 12 through 38. In Figure 12, the distribution of 
the NRCS values for the entire image is presented. The values in the 
image are range from +34 to -36 and have a mean of -15 dB and 
illustrates the diverse of the clutter groups contained within a typical 
airport clutter scene. The histograms shown in Figures 13 through 38 
are plotted as a function of incidence angle, as listed in each figure 
in parenthesis along side the clutter type. In the incidence angle 
range where scattering coefficients are stable, (40° to 60°) 10° 
incidence angle bins were used. In that range where sharp changes in 
scattering coefficient value are expected, (0°) 5° bins were used. The 
incidence angle is equal to 90-depression angle. 

Figures 13 through 17 present the percent of occurrence of NRCS 
values of residential clutter for five incidence angles (90°-depression 
angle) ranges. The mean NRCS value changes slowly with the incidence 
angle in the range from 40° to 75° and fall within a 2.5 dB interval. 
There is also very little change in the shape of the distribution. In 
all cases, the shape of the distribution at the high end is nearly 
identical, while the weak returns are much more sensitive to changes 
with incidence angle. The insensitivity of the strong returns may be 
attributed to man-made targets. The modulation of the weak return may 
come from changing tree/ground interactions. The width of the 
distributions also varies with incidence angle, but the change was not 
systematic. The results for 'forest 1 clutter are presented in Figures 
18 through 23. In general, their cross-sections are 2 to 3 dB lower 
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than those from the residential areas and decreased approximately 4 dB 
with increasing incidence angle (decreasing depression angle). The 
forest clutter distributions are more symmetrical than those of the 
residential clutter, do not change shape with incidence angle, and are 
similar to the distributions obtained using the SAR archival data 
discussed in Report I. 

Unlike the data acquired from the SAR archive, the residential and 
forest areas surrounding the Philadelphia Airport were not similar. The 
distributions of the residential clutter are broader and shorter than 
those of the forest clutter areas and the mean of the residential 
clutter is consistently higher than those of the forest clutter for all 
incidence angle groups. The coefficient of variation of the residential 
areas are generally a factor of ten greater than those of the forested 
areas. This may be due to the diversified nature of the residential 
area. The residential areas not only have trees; they also include 
bright hard targets such as houses and cars, and very low return areas 
such as lawns and streets. 

Figures 24 through 30 present histograms of the scattering 
coefficients of grass for various incidence angle ranges. The mean NRCS 
was almost constant at -24 dB for incidence angles from 40° to 75°. A 3 
dB reduction occurred when angles increased from 75° to 85° region. 

This tendency was typical for distributed targets. Interestingly, grass 
scattering cross-section were only 2 to 3 dB lower than those for 
forest, a difference which is significantly smaller than predicted by 
published data. There is always an uncertainty in the clutter 
characterization when based on topographical and use planning data. 

This, of course, could be corrected with ground truth. Areas which were 
classified as grass may include vegetation which produces much larger 
backscattering cross sections than those typically associated with 
grasslands. Forested areas may also not have been heavily canopied, 
hence would produce a weaker backscatter. The areal coverage of forest 
in the Philadelphia Image was not as extensive as that seen in archival 
imagery, and the character of this forest clutter is probably different. 
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The coefficient of variation for the grass clutter groups is stable 
and the shape of the distribution of values changes very little in the 
40° to 75° incidence angle range. In the 75° to 85° range the shape of 
the distribution is similar to that of the other clutter groups, 
however, these were provided with a much smaller number of pixels, 
because the majority of the pixels had values too close to the minimum 
measurable NRCS (Figures 31 and 32). 

The results of the analysis of the backscatter from water are 
presented in Figures 33 through 38. The source of the water is the 
Delaware River. Inland bodies of water are small and are typically 
sheltered, reducing the surfacing roughening effect possible because of 
winds. Water returns were expected to be low and in a range from -30 dB 
to -45 dB. The minimum measurable NRCS in the image is -36.17 dB. Of 
the six incidence angle ranges examined, 97% to 99% of the water data 
were at or below the minimum value. 

The runway areas of the image produced much the same results as the 
water areas. Experimental backscatter coefficients for asphalt and 
concrete surfaces at 75° to 85° (the location of the Philadelphia 
Airport in the scene) range from -35 dB to -40 dB. Analysis results 
show that all pixels had values of -36.17 dB or less. 

The variation of the magnitude of the scattering coefficients with 
incidence angle are presented in Figures 39 through 42 for each of the 
distributed clutter groups examined in this image. For the three major 
clutter groups observed in this SAR image we would expect to find that 
the backscatter response exhibits a plateau in the middle angle region 
which is then followed by a reduction in backscatter intensity with 
increasing angle. The forest and residential backscatter was 
essentially insensitive to increasing incidence angle in the region from 
45° to 75° while grass started to decrease ground 70°. 

The water backscatter response is constant at -36 dB from 45° to 78° 
because it is at or below the noise floor of the measurement system. 

The sharp rise at 80° incidence angle may be due to system noise 
effects. Most system noise was removed from the image, but it was 
statistically not possible to remove all of the system noise without 
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seriously degrading the image. Any system noise that remained, 
especially at far range, was boosted by the radiometric correction and 
in this case, may have produced scattering coefficients with large 
errors. 

A fourth order polynomial was used to describe the backscatter 
response of forest, residential, and grass with angle. These fits are 
shown in Figures 43 through 45. The polynomial provides a good fit to 
the forest and grass clutter, but a less satisfactory for the 
residential data. At the middle incidence angles (48° to 65° incidence 
angle), the polynomial fits the residential data well, but falls too 
quickly at the larger angles (68° to 78°). An updated antenna pattern 
was used in the radiometric correction of the incidence angle data used 
in these fits. 

A number of man-made targets were analyzed which were located in the 
vicinity of the airport. A chart of these targets and their 
corresponding radar cross-section value is presented in Table 4. The 
RCS (a), in square meters of the target was determined by: 

N 

a = r a r y • Ejt? (Eq. 1) 

where cr? = backscattering coefficient for an individual pixel 
N * number of pixels in extended by a target 
r a = azimuth pixel spacing 
ry = range pixel spacing 

Most of the hard targets analyzed around the Philadelphia Airport have 
RCS values of 20 dBsm or greater with the exceptions of parked 
automobiles which have values from 5 dBsm to 15 dBsm. It is noted that 
even at very large incidence angles (75° and above) many large hard 
targets have RCS values above 35 dBsm. The RCS of docks, increased with 
increasing incidence angle and correlates with the response of a 
dihedral formed between the smooth water and a vertical concrete 
structure. Among the other hard targets no clear trends exist because 
of differences in effective areas and aspect angles of the targets. 
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Distributions of man-made clutter are presented in Figures 46 
through 53. These distributions are of normalized radar cross-section 
values (NRCS), which are calculated by dividing the radar cross-section 
values discussed above by the effective area of the target. Figures 46 
through 50 present the distributions of non-vehicular hard targets. The 
airport terminal building and refinery tank distributions seem to have 
some common shape characteristics. All have steeply increasing trends 
on the low cross-section side of the distribution and a more gradually 
decreasing trends with large cross-section. The building distribution, 
Figure 47, show an additional flatter trend at small NRCS values which 
the other distributions do not have. Similar trends in their 
distributions may be caused by similar geometries. The buildings, 
terminals, and refinery tanks can probably be represented by a dihedral 
formed with a conducting vertical side and dielectric horizontal side. 
All three of these objects probably have some top hat like structures on 
their roofs as well. The distribution plot for the docks, Figure 49, in 
the image has a different shape than the three previous hard-targets. 

The docks produced distributions which are symmetric. The docks 
represent a unique hard-target geometry. They are a flat conductive 
plane sitting up above a conductive surface separated by several 
conductive or dielectric poles. This structure is complicated and not 
well-modeled by a simple dihedral. The distribution for bridge clutter, 
Figure 50, shows characteristics of both types of hard targets. The 
histogram for the bridge clutter appears to be more symmetrical than 
that of the buildings, terminals, and refinery tanks, but the sides of 
the distribution are much steeper than that of the docks. 

The histograms for the vehicular hard targets are presented in 
Figure 51 through 53. The serrated appearance of the plane and barge 
plots is due to the small number of pixels used in the distributions. 

All three cases seem to separate themselves into a group of values at 
the image minimum (in the system noise) and a group of values at a 
larger scattering coefficients. The pixels at the image minimum are 
pixels containing background clutter and weak vehicle scatter. Pixels 
at larger RCS values are those from the dominant vehicle scatters The 
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distribution for the parking lot has three main lobes, one at the image 
minimum, one around -27 dB and one around -4 dB. The sources of the two 
higher lobes in the distribution may be automobiles at different aspect 
angles or a combination of automobiles with other parking lot features 
such as lightpoles. 
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VII. CONCLUSION 


The target-of-opportunity image of the Philadelphia International 
Airport was the first image in the NASA windshear program to be 
collected at small depression angles (7.8°). As such, this image 
represented a unique opportunity to expand the knowledge of the 
backscatter characteristics of clutter to a small depression angles. 

Probably the most distinguishing feature of this SAR image is the 
large areal extent of clutter in which the backscatter was so weak that 
it is probably of no consequence to the ability to detect microburst 
windshear in a ground clutter environment and that the backscatter was 
so weak as to not be measurable (i.e., scattering coefficients were 
smaller than -36 dB) At depression angles smaller than about 13° (which 
represent an area equal to about one half width of the image) virtually 
all distributed targets were buried in system noise. Even at the 
largest depression angles the distributed targets such as forests, 
fields, water, and residential rarely had mean scattering coefficients 
greater than -10 dB. This is consistent with trends observed in the 
archival SAR data which was reported in the Phase I report. 

As found in the Phase I statistical analysis, forest and residential 
clutter are very similar. Residential areas produces backscatter which 
is approximately 2 dB greater than for forest clutter. Both show an 
insensitivity to incidence angle in the range from 45° to 82° and, in 
both cases, the backscatter is moderately weak with scattering 
coefficients of about -20 dB. Overall, only a very small percentage of 
the image is composed of high returns, about 3% of the returns were 
greater than -5 dB. Sources which produced the largest cross-sections 
were largely confined to the airport grounds and areas highly 
industrialized. Unlike the original archival data, where more of the 
high backscatter returns are from areas outside of the airport grounds, 
this scene more closely emulates an airport in an urban setting and a 
more complicated clutter environment, which is a good case to be 
exercised in the NASA Radar Simulation Program. 
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TABLE 2. IMAGE COMPOSITION AREAL ANALYSIS 
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TABLE 3. STATISTICAL SUMMARY OF 

BACKSCATTER CROSS-SECTIONS FOR 
VARIOUS CLUTTER TYPES NEAR THE 
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Refinery Tanks -4.98 -36.17 18.73 0.2894 x 10-1 10.05 


Table 4. Hard Targets Represented as o 


Identifier 

Region 

0 

a 

Effective Area 




(dBsm) 

(m 2 ) 

HI 

Dock 

60.37 

33.84 

7573.82 

H2 

Dock 

63.49 

35.14 

7589.38 

H3 

Dock 

65.20 

36.15 

7690.46 

H4 

Dock 

72.26 

38.74 

7612.70 

H5 

Dock 

72.88 

38.39 

7636.03 

H6 

Dock 

73.55 

34.35 

7394.98 

H7 

Dock 

74.25 

41.08 

7737.12 

H8 

Bridge 

76.42 

27.28 

777.60 

H9 

Bridge 

76.72 

33.30 

777.60 

H10 

Bridge 

76.87 

37.15 

777.60 

H1 1 

Bridge 

77.10 

32.77 

777.60 

H12 

Football Stands 

71.05 

26.19 

4611.17 

H13 

Arena 

70.19 

38.89 

42698.02 

H14 

Bridge 

75.73 

33.71 

777.60 

H15 

Bridge 

78.79 

20.06 

373.25 

H16 

Bridge 

78.92 

24.85 

334.37 

H17 

Terminal 

81.40 

37.80 

1407.46 

H18 

Terminal 

81.75 

43.59 

2099.52 

H19 

Plane 

81.17 

24.84 

178.85 

H20 

Plane 

81.84 

22.08 

124.42 

H21 

Warehouse 

80.04 

34.70 

5116.61 

H22 

Barge 

78.17 

25.32 

326.59 

H23 

Barge 

78.75 

31.15 

2122.85 

H24 

Refinery Tank 

53.63 

29.08 

6772.90 

H25 

Refinery Tank 

55.56 

33.08 

4712.26 

H26 

Refinery Tank 

56.88 

32.99 

4758.91 

H27 

Refinery Tank 

58.03 

32.85 

4813.34 

H28 

Hangar 

80.28 

36.44 

816.48 

H29 

Hangar 

80.91 

39.17 

2177.28 

H30 

Ship 

80.91 

37.83 

1458.22 

H31 

Building 

69.46 

30.92 

11749.54 

H32 

Building 

78.84 

42.32 

19813.25 

H33 

Building 

82.44 

42.49 

4922.21 

H34 

Barge 

82.25 

31.56 

373.25 

H35 

Barge 

82.32 

36.30 

349.92 

H36 

Building 

81.34 

47.24 

3654.72 

H37 

Parked Auto 

68.14 

8.19 

23.33 

H38 

Parked Auto 

68.25 

5.62 

23.33 

H39 

Parked Auto 

68.34 

-17.17 

7.78 

H40 

Parked Auto 

68.43 

9.62 

23.33 

H42 

Parked Auto 

68.47 

14.65 

23.33 
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Parameter 

Value 

Altitude 

Depression Angle 
Distance to Touchdown 
Point 
Resolution 
Pixel Spacing 
Image Size 
Polarizations 
Frequency 

1524 m 
8° 

10844 m 

2.7 m * 3.0 m 
1.62 m* 2.4 m 
1 2 km * 9.8 m 
VV, VH, HH, HV 
X-Band 



Figure 1 . SAR Imaging Parameter 


28 







29 








+ 


Near Range 


Far Range 


Figure 3c. Quad-Swath Mode 


Figure 3. P-3 SAR Swath Modes 
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Figure 4. SAR Data Collected of Philadelphia Airport and 
Surrounding Area on 7 October 1988 
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Bar Chart Presentation 
of Means and Standard Deviations 
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Figure 5. Bar Chart Presentation of Scattering Coefficients for a 
Variety of Clutter Types Near the Philadelphia Airport 
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Figure 6. Histogram of Radar Scattering Coefficients for the 
Philadelphia Image 
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Figure 7. Philadelphia SAR Image Thresholded at -30 dB 
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Figure 8. Philadelphia SAR Image Threshholded at -20 dB 
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Figure 9. Philadelphia SAR Image Threshholded at -1 0 dB 
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Figure 10. Philadelphia SAR Image Threshholded at 0 dB 
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Figure 1 1 . Philadelphia SAR Image Threshholded at 1 0 dB 
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Figure 12. Distribution of Radar Scattering Cross Sections in the 
Philadelphia SAR Image 
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Figure 14. 
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Figure 17. 
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Figure 19. 
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Figure 25. 
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Figure 40. incidence Angle vs Scattering Coefficient Plot for 
Residential Clutter 
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Figure 41 . 
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Figure 43. Polynomial Fit to Forest Data 


70 


Scattering Coefficient (dB) 

- 50.0 - 40.0 - 30.0 - 20.0 - 10.0 


Residential Clutter 


O 

O 



j 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 t i | r 

40.0 50.0 60.0 70.0 80.0 

Incidence Angle (degrees) 


Figure 44. Polynomial Fit to Residential Data 
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Figure 45. Polynomial Fit to Grass Data 
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Abstract 

A new triband polarlmetrlc synthetic aperture radar 
(SAR) has been developed for remote sensing applica- 
tions and has been Installed In a U.S. Navy P-3 
aircraft. Pulses are transmitted at either X-, C-, or 
l- band. They may be transmitted at either horizontal 
(H) or vertical (V) polarization and received at either 
horizontal or vertical polarization. Transmit and 
receive polarization, as well as frequency band, may be 
switched pul se-to-pulse. Up to four different 
waveforms may be Interleaved, e.g., X^h Lhv *VV ^VH or 
*1111 X|IV *VH *VV (ful 1-polarlmetric mode). A digital 
synthesizer produces the FM-chlrp pulse. The pulse Is 
emitted from a new triband antenna that Is kept perpen- 
dicular to the aircraft ground track by means of a 
three-axis, servo-driven positioner which compensates 
for aircraft pitch, roll and yaw. Received pulses are 
digitized, preprocessed and recorded to high-density 
digital tape. Test flights began In late October, 19B7 
and preliminary Imagery has been obtained. 


Introduction 

A SAR system, called the P- 3/S AR f has been 
Installed In a U.S. Navy P-3 aircraft and operates at 
X - , C- and L-bands with the capability of recording HH, 
HV , VH and VV polarizations. An earlier description 
was provided [1], The system Is a fully focused side- 
looking SAR capable of looking out either side of the 
aircraft. Operating center frequencies are 9.35 GHz, 
5.30 GHz and 1.25 GHz. The system data are recorded in 
several forms. Hlgh-denstty digital tape (HOOT) Is the 
primary recording medium. The data are also recorded 
on photographic film for subsequent optical processing. 
A real-time Image formation processor records Image 
data on heat-developed film or paper. 

Figure 1 Is a photograph of the P-3 aircraft. The 
radome aft of the wing contains the antenna and 
pedestal . 



Op erating M odes 

The P-3/SAR has five operating modes, summarized in 
Figure 2. For each mode, the figure illustrates the 
number of frequcncy/polarlzatlon combinations, the 
pulse repetition frequency (PRf) per comb! fiat Ion, the 
number of range bins per pulse, the presum factor (data 
are presummed via a presum filter before being written 
to tape) and the corresponding data recording rate, 
tape speed and time to record a full tape. The rate of 
recording data to tape Is 

r = 2nf R Nb/P 

where 

n = number of f requency/pol arl zat Ion combinations 
f R = PRF 

N = range bins per pulse 
b = bits In the A/O converter = 6 
P “ presum factor. 

The factor of 2 Is to account for both I and Q. 

The Single-Swath Multiplex Mode employs any four 
frequency/polarization combinations, except that (1) at 
least two different frequency bands must be used 
because of the transmitter duty- cycle limitation, and 
(2) each band must have a uniform PRF within the over- 
all waveform. The PRF per combination Is 2 kHz at a 
ground speed of 350 knots. (PRF Is proportional to 
ground speed producing a constant value of pulses per 
unit ground distance). Four thousand ninety-six range 
bins per FM-chlrp pulse are processed. These are pre- 
summed to a data rate one-sixth that which would result 
from no presuminlng. An example for this mode Is shown 
In Figure 3. 

The Double-Swath Multiplex Mode Is similar except 
that for each pulse, the number of range bins Is 
doubled to B192. The number of f requency/polarl zatlon 
combinations Is halved, to two, and the presum factor 
remains the same at six. Thus, the data recording rate 
Is the same. An example Is given In Figure 3. 1 he 

Quadruple-Swath Multiplex Mode Is completely analogous 
with one frequency/polarization combination, 16,381 
range bins per pulse, a presum factor of six and the 
same data recording rate. An example Is shown In 
Figure 4. 

The Single-Swath Polarlmetrlc Mode Is similar to 
the Single-Swath Multiplex Mode exenpt that all pulses 
are from the same frequency hand and nil four polariza- 
tion combinations (Mil, UV, VII, VV) are employed. IMs 
enables the full polarization scattering matrix to be 
recorded. Because of the transmitter duty-cycle limit 
of two percent , the PRF per combination Is limited to I 
kHz instead of 2 kHz. Hie presum factor Is reduced to 
three, again resulting In the same data recording rate. 
An example Is given In Figure 4. 


FIGURE 1. P 3 AIRCRAFT 
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MODE 
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- n 

PRO* per 
Frcr^Pol 
Combi nation 
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( 117 ) 
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Recoicling 
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2irl|| NO 

" P 

(MtiH-J-n-c) 

Nominal b 
lope Speed 
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tape time 
1- lAr 

Single Swath 
Multiplex 

4 

2000 

4096 

6 

05.5 

(71) 

ct C 
|8G 2) 

35 4 
(25 3) 

Double-Swath 

Mulliplex 

2 

2000 

8192 

6 

05.5 

(71) 

61 6 
{062) 

35 4 
(25 3) 

Quadruple Swath 
Multiplex 

1 

2000 

16384 

6 

05.5 

(71) 
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(OG 2} 
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(25 3) 

Single Swath 
Polarimeltic 

4 C 

1000° 

4096 

3 

055 

(71) 
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(8G2) 

35 4 
(25 3) 
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No Presum 

1 

2000 

4096 

1 

903 
(10C. 5) 

92 4 
(1293) 

23 6 
(169) 


a - PRF is proportional lo ground speed. Values given are appropriate lor 350 knots 

b - Actual values, given in parcnlheses, are slightly dilterenl due to "second order" ellecls such as record lotmnl. 


auxiliary dala, etc. 

C- Transmitter duty cycle is 4000 pulses/sec x 4 psec/pulse - 1.6 percent. Limil is 2 percent, preventing higher Pill' 


b - bits in A/D - 6 
n t - no. tracks - 28 
p - tape density - 36,000 bits/in Track 
I - tape Icnglh - 10.888 It - 130.65G in 


FIGURE 2. OPEItAIING MODES 
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FIGURE 3. MODE TIMING DIAGRAM 


FIGURE 4. MODE TIMING DIAGRAM 
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The Single-Swath No-Presum Mode provides a capabil- 
ity for recording all of the Doppler data spatially 
filtered by the antenna beam. No additional azimuth 
filtering Is applied to this data In the preprocessor. 
In this case, only one frequency/polarlzatlon combina- 
tion can be used and the presum factor Is one. The 
data recording rate Is 1.5 times that for the other 
modes. 


Subsystem Configuration 

Figure 5 shows a simplified block diagram of t fie P- 
3/SAR system. 



mio o 

H> NAOC 


FIGURE 5. SYSTEM BLOCK DIAGRAM 


The Inertial Navigation System (INS), a Litton UN- 
51, performs several important functions for the radar 
system. It steers the aircraft via the autopilot 
during mapping, It points the antenna to 90 degrees 
from the ground track, It provides aircraft rotation 
data to the antenna glmbal systems for compensation of 
pitch and roll, and It provides aircraft lateral trans- 
lation sensing for motion compensation of the radar 
data. Integrated INS three-axis accelerations are 
digitized In the INS Interface and are input to the 
data concentrator. This computer resolves these Inputs 
from the INS coordinate system into the radar coordin- 
ate system (along-track , cross-track, vertical) and 
then computes the velocity component along the radar 
beam. This beam velocity Is scaled to radar wavelength 
and converted to phase corrections. 

The linear FM chirp waveform Is generated within 
the Recel ver/Excl ter/STALO using a digitally synthe- 
sized waveform generation technique described by 
Callaway , et al. In another paper In these Proceedings. 
The chirp data points are stored In a PROM, then loaded 
Into a RAM. The data points are read out of the RAM at 
the clock rate to produce the digital chirp waveform. 
This waveform is converted to an analog signal, then 
filtered and up-converted to the desired RF carrier 
frequencies for transmission. One feature of this 
technique is the capability to alter the chirp waveform 
data points In the PROM to compensate for range 
dimension hardware phase errors. 


The chirp waveform Is Input to the transmitter 
where It is amplified by a traveling wave tube (1WT) 
amplifier. The X-, C- and l -band IWIs are manufactured 
by Hughes, Litton and Varlan, respectively. A four 
microsecond pulse Is transmitted for each band. The 
pulse repetition frequency (PRF) Is proportional to 
ground speed to provide constant azimuth spatial 
frequency sampling. PRF Is 8 kHz for a ground speed of 
350 knots. 1 lie maximum transmitter PRT for any band Is 
4 kHz with the maximum transmit duty cycle equal to 2 
percent for each band. 

Ihe amplified chirp waveform Is passed to the 
Transml tter/Duplexer where a hlgh-power switch selects 
the transmit polarization. 

A triband (X-, C-, L-bands) dual -pol ar! zat Ion 
antenna, developed by Chu Associates of Littleton, 
Massachusetts, Is used for this system. 1 he X- and C- 
band capability is provided by a dual, feed-on- focus 
shaped reflector based on the work of L.J. Chu (2), 

The L - band capability Is provided by a crossed-clemcnt 
Yagl array mounted on top of the X-C band reflector. 
This antenna is mounted on a new three-axis glmbal 
pedestal developed by Klntec Corporation of St. 
Petersburg, Florida, This glmbal system provides 
stabilization of the antenna with attitude reference 
data from the INS. 

8ackscattered signals are collected by the antenna 
and passed through the Transmltter/Uuplexer where a 
switch selects the receive polarization. Ihe 
Transml tter/Duplexcr Includes a transmitter power mon- 
itor capability which continually monitors average 
transml tter power. In addition, an Injection port for 
amplitude-controlled synthetic targets Is provided. 
This gives an end-to-end test capability lor pre-flight 
and post-flight testing and system calibration. 

The signal Is then sent to the receiver which pro- 
vides down-conversion of the X- and C-band RF channels 
to the l -band IF, band-limits the signals to support 
the desired range resolution and then amplifies the 
signals with minimal addition of noise. A digital 
phase shifter and frequency synthesizer within the 
receiver provide the motion compensation correction for 
the received signals based on the Inputs from the INS. 

The Recel ver/Excl ter/S l AID provides complex I (In- 
phasc) and Q (quadrature) analog video signals to the 
Digital Preprocessor, and real analog-video signals to 
the Cttl/fl hn recorder and clutterlock units. 

The Digital Preprocessor performs A/D conversion of 
the I and Q analog video signals, buffers these signals 
to reduce the data rate, performs 6:1, 3:1 or no pre- 
suming of the Input signals depending on the radar 
mode, and then formats the data for high-density digi- 
tal tape recording (MODI R ) on the Honeywell AM/USM-XX 
28 track tape recorder. The two-channel A/D converter 
required to digitize the complex radar video utilizes 
two 6-l>lt high-performance A/D converter circuits built 
at ER1M. 

The AN/USH-XX HDOTR 28 tracks are utilized as fol- 
lows: 24 tracks for four 6-blt complex video words, 
one track for auxiliary data (128 16- bit words per 
frame), one track for frame synchront zat Ion, and two 
tracks for error correction. Four thousand ninety-six 
complex words make up a frame of data. 

The Digital Preprocessor also provides p re summed I 
and Q digital video to the Rcal-llme Processor (RIP). 
The RIP performs single-channel, Image-formation pro- 
cessing In real time. Ihe output Is displayed on a 
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VISOR hardcopy thermal paper recorder, 
vldes capability for real-time check 
navigation and also a coarse check 
operation. 


This RTF pro- 
of the general 
on the system 


The CRT /Film Recorder provides a capability for 
recording one channel of video data on silver halide 
film for subsequent ground processing in an Optical 
Processor. 


The clutterlock unit uses a low- f rcquoncy sprclintn 
analyzer to measure the Doppler frequency centroid of 
the radar returns spatially filtered by the antenna 
beams. A cross-track velocity analog Is developed from 
this Doppler frequency centroid; It Is input to the 
Data Concentrator velocity filter to remove INS 
velocl ty biases. 

The system computer consists of two HP 9000 Series 
3?0 computers and peripherals to (1) provide operator 
control of the radar system signals, (2) perform the 
motion compensation computations, (3) provide an 
operator/display Interface, and (4) concentrate data 
for auxiliary data recording and other functions. The 
Control Computer and associated monitor serve as the 
operator interface. This computer does all of the 
slower functions such as mode control, line navigation, 
system monitoring, etc. The Data Concentrator performs 
the motion compensation computations, converts the 
roll, pitch and yaw INS measurements into antenna stab- 
ilization and pointing commands, filters and combines 
the different sources of velocity measurements (INS, 
GPS and clutterlock), and assembles the auxiliary data 
for recording on the HDDTR. 

The Time Code Receiver uses the GOES satellite data 
as a time base for comparing radar events with other 
data collection equipment events. This is distributed 
to other systems within the P-3. 

The Global Positioning System (GPS) receiver pro- 
vides an independent measure of aircraft position and 
velocity. The aircraft position (latitude and longi- 
tude) Is used for manual update of the INS to Improve 
the general navigation capabilities. Ilie GPS veloc- 
ities are resolved Into along-track and cross-track 
velocities. The along-track velocity is used for set- 
ting the system PRF. The cross-track velocity is used 
as an Input to the velocity filter to reduce INS 
biases. These Inputs will be used when GPS data Is 
available. 


The system output products are (1) high-density 
digital tapes from the Honeywell AN/USH-XX HUDIR, (2) 
signal film cassettes from the CRT/Film Recorders, (3) 
a Control Computer printout which Includes mode param- 
eters, pass parameters, flight geometry and operator 
comments, and (4) 9-1nch paper format real-time Imagery 
from the VISOR recorder. The high-density digital 
tapes are played back In an ERIM facility where 
computer compatible tapes are generated for subsequent 
image-formation processing. The signal film data Is 
chemically processed, then used In the optical proces- 
sor to expose an Image film. 

Images are processed In one of three ways: 

Optical -Optical : Raw video data from the on-board 
optical film recorder are processed via Fourier optics 
to produce Imaging. This Is the fastest processing 
method. 


record, and then processed via Tourler optics to pro- 
duce Imagery. This can provide a quick look at the 
contents of an entire HOUf. 

yia.lta!-PJgUaLs Oata from the HUIM aie processed 
via computer to produce Imagery. jlils Is the most 
precise method for analyzing a specific frame but Is 
not generally suited to an overview of an entire 
f I Ight. 


Resu l ts 


The first test flight of the P-3/SAR flight, 
designated as an engineering flight, occurred ott 28 
October 1987. Ihe aircraft took off from Willow Run 
Airport near Ann Arbor, Michigan and was flown to the 
Straits of Mackinac and back twice, at an altitude of 
15,000 feet. At the time, permission to transmit at L- 
and C-bands had not yet been received. Iherefore, 
transmission was planned at X-band only. Data were 
collected at Xgg and X^v. 

Tlgure 6 shows "digital -dlgl tal “ Images of Saginaw, 
Michigan, taken at and X^y. For each Image, the 
gray scale was chosen to reveal maximum detail; thus, 
one cannot compare relative brightness between the two 
Images. In general, the Xpg Image contains many more 
sharp, specular reflections, for example In the region 
of the Saginaw Fair Grounds. It Is Interesting tn note 
the return from the railroad tracks at both X H u and 
Xpy. Hie long north-south portion of the tracks Is 
essentially parallel to the aircraft flight path; there 
Is also a set of tracks approximately perpendicular to 
the flight path. At Xp||, where the transmitted and 
received electric vectors are parallel to the flight 
path, the "parallel" tracks produce a brighter return 
than the "perpendicular” tracks. At Xgy, the perpen- 
dicular tracks are brighter than the parallel tracks. 
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Digita l-Optica l : Data from the blgh-denslty data 
tape (HUD!) are transformed to video data on film, In 
the same format as that produced by the optical film 
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APPENDIX B 
SAR CALIBRATION 

In this Appendix the calibration of the Philadelphia image is 
discussed in detail. When a SAR system is operated in the linear region 
there exists a one-to-one correspondence between the image intensities 
of targets in a scene and the radar cross-sections of those targets. To 
account for radar system gain effects and the benefits gained in 
optimizing the recording process, the data undergoes scaling as it is 
collected. One purpose of absolute calibration is to determine this 
constant and to rescale the image intensities values so they now 
represent absolute scattering cross-sections. If the cross-sections are 
normalized by the pixel area than they are referred to as normalized 
radar cross-sections (NRCS) or scattering coefficients. 

Since there were no reflector arrays in the Philadelphia image, we 
used a calibration scene collected at the Patuxent Naval Air Test 
Center, Maryland, two weeks prior to the Philadelphia image to calibrate 
the data. This can be done because the P-3 SAR system transfer function 
is continuously monitored throughout the data collection. 

As discussed in the first SAR data report, the linear relationship 
between a, the expected radar cross-section of a reference target, and 
PI, the measured image intensity from a radiometrical ly corrected image 
of this target is given by: 


Pr- P N Eq. (B-l) 

a - 

K 

where P^ is the system noise and K is the system gain constant. K is 
determined from the intercept of the calibration function with the Pj 
axis. 

The radar scattering cross-sections for each of the triangular 
trihedral reflectors used in the calibration array may be accurately 
calculated using the equation provided in Ruck et. al., (1976) 
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a{6 ,<t>) = a cos0+sin0(costf+sin$) - ^ 

o 

X [cos0+sin0(cos0+sin0)] 

Eq. (B-2) 

where 6 is the angle between the vertical crease of the trihedral and 
the direction of illumination, <p is the offside angle to the trihedral, 
a is the length of the trihedral along a crease, and X is the wavelength 
of the radar. In this calibration array, the trihedrals were placed so 
that they were flat on the ground and their aperture faced the look 
direction of the radar. In this case, the angle 6 is simply equal to 
the incidence angle to any reference trihedral is given as: 

6 = cos" 1 _!!_ Eq. (B-3) 

Rs 

where H is the altitude of the radar and Rs is the slant range to the 
reflector. The slant range to the reflector is given by: 

Rs = (Td * 150.00) + (Prl -1) * Rps Eq. (B-4) 

where the quantity Td is the time delay to ground cells in the near 
range of the image, Prl the peak record location, and Rps the pixel 
spacing. The values of a, crease length, were 0.3 m, 0.45 m, 0.6 m, 

0.75 m, 0.9 m, and 1.2, for the reflectors used in the calibration 
analysi s . 

The image intensities of each of the calibration reflectors were 
measured by positioning a 32 x 32 pixel window around the apparent peak 
of each reflector. Upsampling was performed to improve the accuracy of 
the estimate of the energy associated with the reflector. The 
upsampling was accomplished by taking the Fourier transform of the data 
contained in the 32 * 32 pixel window. In the frequency domain the 
bandwidth of the data is determined by the sample size of the image 
data. By 'zero-padding' the data in the Fourier domain the information 
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content was preserved but pixel spacing in the image domain was reduced. 
Given the upsampled data, the total reflector energy was calculated. 
These values were then radiometrical ly corrected and the effects of 
local system noise and background clutter were removed. The average 
intensity of the sum of the clutter and system noise was calculated by 
finding the mean intensity value of a small area near, but not 
interfered by, the target of interest. The amount of background energy 
which contributes to the reflector energy, P B , is given by: 

P B = ((V B )**2)*r a *r r 

where V B is the average background amplitude near the reflector, r a is 
the azimuth resolution, and r r is the range resolution. Subtracting 
this from the total 3 dB image intensity values removes the background 
clutter and system noise, Pn. 

Next, the effect of antenna gain pattern and range fall -off were 
removed from each image intensity value. In essence, we performed a 
radiometric correction to each reflector energy identical to that which 
is performed on the entire image. The transmit and receive antenna 
patterns act as weighting factors in the image. Each pixel of the image 
is multiplied by a scale factor which is the product of the transmit and 
receive antenna pattern. The pattern used to correct the Philadelphia 
Airport image is illustrated in Figure B-l. To remove the effects of 
antenna pattern one needs only to multiply the measured image intensity 
value by the inverse of the transmit and receive antenna patterns. 

Range fall -off is removed by multiplying the measured intensity values 
by the square of the slant range to the point of interest. An R 2 factor 
is employed instead of an R^ because the SAR processor removes two 
orders of the range fall -off effect as the data is convolved with its 
match filter in azimuth. The data is then normalized to a reference 
range, which is set to be slant range at the antenna boresight. The 
corrected reflector energy Prc is then: 
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(P r - P B ) * R s **2 


(R ref **2)*G t *G r 

where P r is the total 3 dB energy, Pg Is the total background energy, R s 
is the slant range to the reflector, R re f is the reference slant range, 
G t is the transmit antenna pattern gain, and G r is the receive antenna 
pattern gain. 

The normalized reflector image intensity P rc is then plotted against 
the expected RCS of the reflectors and the system gain constant is 
determined from the intercept of a least squares fit line through the 
data with a slope of one. In the case of the Philadelphia Airport 
image, since the calibration was performed on another image, we also 
adjusted for differences in transmitted power, stable local oscillator 
attenuation, processing scale, depression angle, and altitude. The 
calibration data are plotted in Figure B-2. 
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APPENDIX C 

STATISTICAL FORMULATION 


In this appendix, we describe the calculation of the statistics used 
to describe the scattering properties of the various clutter types. The 
mean p of a data set is given by: 



N 

E x 


Eq. (C-l) 


where N is the total number of data values and x^ ( represents an 
individual data values. 

The variance a 2 is given by: 


a 2 = 



E (xj - p) 

i=l 


Eq. (C-2) 


The standard deviation is defined as the positive square root of the 
variance. The coefficient of variation is given by the ratio of the sum 
of the standard deviation plus the mean to the mean of the data set. 


r o + V- 

L v 


Eq. (C-3) 
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APPENDIX D 


DETAILS OF CLUTTER ANALYSIS SURVEYS 

In Figure D-l, the coded locations of sites used in the statistical 
characterization of clutter in the Philadelphia SAR image are provided. 
The letter indicates a clutter type; such as, I for industrial, W for 
water, H for hard targets, G for grass, M for man-made cultural sites, R 
for runway, U for urban, F for forest, and Re for residential. In Table 
D-l, the results of the statistical analysis of these individual clutter 
sites are provided. In Table D-2, the results for general clutter 
categories are provided. The clutter type and range of incidence angles 
are included in the file names. These data were generated from those 
described in Table 2. 
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Figure D-1 . Locations Where Clutter Analysis Surveys of the 
Philadelphia Airport SAR Scene Conducted 
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Table D-1 . Statistical Analysis Results for Individual Clutter Groups 
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